ANALYTICITY AND UNIFORM STABILITY OF THE INVERSE 
SINGULAR STURM-LIOUVILLE SPECTRAL PROBLEM 



ROSTYSLAV O. HRYNIV 

Abstract. We prove that the potential of a Sturm-Liouville operator depends ana- 
lytically and Lipschitz continuously on the spectral data (two spectra or one spectrum 
^vq and the corresponding norming constants). We treat the class of operators with real- 

^ ' valued distributional potentials in (0, 1), s € [0, 1]. 

00 

^| 1. Introduction 

p | 1 In this paper we shall establish analyticity and uniform stability of solutions of 

two inverse spectral problems for a certain class of Sturm-Liouville operators on the 
! interval [0, 1]. The (direct) spectral problems to be considered are 

(1.1) -y"(x) + q(x)y(x) = Xy(x), 

subject to some boundary conditions, where A G C is the spectral parameter and q 
is a real-valued potential that might be regular (i.e. integrable) or singular (e.g. a 
distribution). For simplicity, we shall restrict ourselves to the cases of the Dirichlet 
boundary conditions 

(1-2) „(0) = y(l) = 

and the Dirichlet-Neumann boundary conditions 
(1.3) y(0) = y'(l) = 0, 

although other boundary conditions can be treated similarly. (We note here that the 
derivative in ( II. 3p must be replaced with a quasi- derivative if q is singular, see Section [2] 
for details.) We shall always denote by Ai < A2 < • • . the eigenvalues of problem (11. ip . 
ffl~2]) and by fn < fj, 2 < ■ ■ ■ those of problem (fTTTj) . (TOD . 

In 1946, Borg [7J proved that the spectrum of problem (11. ip corresponding to one 
fixed set of boundary conditions, e.g., the Dirichlet (II. 2p or the Dirichlet-Neumann (II. 3p 
ones, does not determine the potential q uniquely. (An exceptional situation where only 
one spectrum — namely, that for the Neumann boundary conditions — determines the 
problem was pointed out in 1929 by Ambartsumyan [3].) However, two such spectra 
already suffice to reconstruct the potential q unambiguously, as follows from the in- 
verse spectral theory for Sturm-Liouville and Schrodinger operators that emerged in 
the early 1950- ies from the work of Gelfand and Levitan [T4], Marchenko [37], and 
Krein [33] and has been extensively developed in many directions since then. This the- 
ory gives an efficient algorithm reconstructing the potential q from the spectra (A n ) n6 N 
and (/U„) ne N and also describes completely the set of possible spectra. For instance, a 
typical result [381 Thm. 3.4.1] reads 
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Theorem A. For real numbers Ai < A2 < . . . and /ii < fi2 < ■ ■ ■ to give all the 
Dirichlet and Dirichlet-Neumann eigenvalues of the Sturm-Liouville problem (11.11) with 
a real-valued q G L 2 (0, 1), it is necessary and sufficient that these numbers 

(i) interlace, i.e., fi n < X n < X n+ i for all n G N, and 

(ii) have the representation 

X n = 7r 2 n 2 + A + a n , fi n = n 2 (n - \) 2 + A + b n , 
where A G E and (a„) ne N and (6 n )neN ore some ^-sequences. 

The induced mapping from the spectral data ((A n ), (/!„)) into the potentials q 
provides a solution to the inverse spectral problem and has been extensively stud- 
ied in the literature. In particular, this mapping is shown to be locally continuous 
in a certain sense, which yields local stability of the inverse spectral problem, see, 
e.g., 0Cl[I2linil2nraE^I32ISlSSlllSMIS3 and the references therein. Here we 
introduce a metric on the set of the spectral data ((A n ), (fi n )) by e.g. identifying such 
data with the triplets [A, (a n ), (&„))g R x f 2 x 4 in the representation of item (ii) 
above. Typically, this local stability states that, for a fixed M > 0, there are positive e 
and L with the following property: if potentials q\ and q 2 are such that ||<Zi||l p (o,i) — M 
and ll^llLp^i) < M and the corresponding spectral data v\ : = ((Ai >n ), (}J-i,n)) and 
v 2 := ((A 2 ,n), (j"2,n)) satisfy ||^i - u 2 \\ < e, then 

(!-4) ||?i - 9211^(0,1) < ^ll^i ~ ^2 1 1 

for a suitable p6 [1, 00]. For instance, local stability results with p = 2 were established 
in pf3jloT)] in the regular case q G L 2 (0, 1), and in for impedance Sturm-Liouville 

operators. The cases p > 2 and p = 00 were treated in [23] and in [T7JI11] respectively; 
earlier Hochstadt in [2T|[22] proved stability if only finitely many eigenvalues in one 
spectrum are changed. The papers [201S8] studied to what extent only finitely many 
eigenvalues in one or both spectra determine the potential, and the latter problem in the 
non-self-adjoint setting was recently discussed in [3D]. In [36], stability of reconstruction 
of general first-order systems from various given data was investigated. Also, stability 
of the inverse spectral problems on semi- axis was proved in [39,49], and the inverse 
scattering and inverse resonance problems on the line and half-line were studied in [TPj 
HP] and p2HHlg2] respectively. 

However, the above results cannot be considered satisfactory, as they refer to the 
norm of the potential q to be recovered and thus specify neither the allowed noise 
level e nor the Lipschitz constant L. Therefore we need a global stability result that 
asserts (11.41) whenever the spectral data V\ and u 2 run through bounded sets ,JV and 
with L only depending on jV . 

Such global stability of the inverse problem of reconstructing the potential q of 
a Sturm-Liouville equation from its Dirichlet and Dirichlet-Neumann spectra was 
recently established by Savchuk and Shkalikov in the paper [ST]. In fact, the au- 
thors considered therein the class of problems with distributional q in the Sobolev 
space 14 7 2 _1 (0,1), s G (0,1]. Two typical and most important examples of singular 
potentials belonging to W 2 ~ (0, 1) only if s < \ are the Dirac delta-function 6(- — a) 
and the Coulomb-type interaction l/(- — a), a G (0,1). Recently, singular Sturm- 
Liouville operators of various types including e.g. operators with distributional poten- 
tials in H^2 _1 (0, 1) [25|I28|I29|I5TH55] and operators in impedance form [UEJIEJIIT] have 
attracted considerable interest, which is partly motivated by their importance for many 
applied problems in classical and quantum mechanics, scattering theory etc. Both the 
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direct and the inverse spectral theory of such singular operators were developed in 
detail; we refer the reader to e.g. [2"5 |l55ll5it l5B] . 

The purpose of the present paper is two-fold. Firstly, we give an alternative proof 
of the global stability of the inverse spectral problem for Sturm-Liouville operators 
with distributional potentials in M^ 2 S_1 (0, 1) for s G [0,1], i.e., including the extreme 
case s = (Theorem 12. ip . In fact, we prove that the inverse spectral mapping is 
analytic and locally Lipschitz continuous. Secondly, we prove analogous results for 
the inverse spectral problem of reconstructing the potential q from the corresponding 
Dirichlet spectrum (A n ) and norming constants (a n ) defined in Sectionl2l(Theorem l2.2l) ; 
the local stability in this setting was established e.g. by McLaughlin in [13]. Our 
approach is completely different from that of the paper [13J (which follows the method 
developed by Poschel and Trubowitz in the book [IB]) and of the paper [57J (which uses 
the modified Priifer angle) and is based on a generalization of the classical Gelfand- 
Levitan-Marchenko method as developed e.g. in [25] . The proof essentially relies on 
the fact that certain sets of sines and cosines form Riesz bases of 1/2(0, 1) with uniformly 
bounded upper and lower bounds [21]. As in [57] we prove first the required results for 
the extreme cases s = and s — 1 and then use the Tartar nonlinear interpolation [58] 
to cover the intermediate s. 

The paper is organised as follows. In the next section we formulate the main re- 
sults and define suitable topologies on the sets of spectral data. Section [3J describes 
the reconstruction method based on the Gelfand-Levitan-Marchenko integral equa- 
tion. Theorem 12.21 on analyticity and local Lipschitz continuity of the inverse spectral 
mapping using the Dirichlet spectrum and norming constants is proved in Section HJ 
and Section [5] studies dependence of the norming constants on two spectra. Finally, 
Appendices lAj El and [C] contain some auxiliary results on Riesz bases, Sobolev spaces, 
and Fourier transforms therein. 

Notations. Throughout the paper, ^fz shall denote the principal branch of the square 
root that takes positive values for z > 0. For a Hilbert space H, we denote by 38(H) 
the algebra of bounded linear operators acting in H. 

2. Preliminaries and main results 

In this section we define explicitly the class of Sturm-Liouville operators to be studied 
and state the main results for the inverse problems. 

2.1. The operators. Given a real-valued distribution q in W 2 ~ 1 (0, 1), we define Sturm- 
Liouville operators T corresponding to the differential expression 

d 2 

dx z 

by means of regularisation by quasi-derivatives suggested by Shkalikov and Savchuk [521 
[53]. We take a real- valued distributional primitive a G L 2 (0, 1) of q and set 

h(y) ■= ~(y' - <ry)' - °y' 

for y G W^O, 1) such that the quasi- derivative := y> _ y j s absolutely continuous 
and l a (y) is in L 2 (0, 1). We then define the operators T D = T D (<r) and T N = T N (o~) as the 
restrictions of l a onto the functions satisfying the boundary conditions y(0) = y(l) = 
and y(0) = y^(l) = respectively. 

Since l a (y) = —y" + qy in the sense of distributions, these operators coincide with 
the classical Sturm-Liouville operators if q G Li(0, 1). Moreover, the operators Tb(cr) 
and T N (cr) depend continuously on cr G L 2 (0, 1) in the uniform resolvent sense [521155] . 
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Therefore Td(ct) and Tn(<t) are the most natural Sturm-Liouville operators associated 
with differential expression dHJ for q G W 2 \0, 1). 

We observe that although the differential expression l a is independent of the partic- 
ular choice of the primitive a of q, the boundary conditions for the operator T N (cr) and 
the norming constants for Td(ct) introduced below do depend on a rather than on q. 
Therefore it is the function a, and not q, that has to be reconstructed in the inverse 
spectral problem for singular Sturm-Liouville operators under consideration. We shall 
call a the regularized potential of the Sturm-Liouville operators Tb(cr) and Tn(ct). 

2.2. Spectral data. It is known [521153] that for a real- valued o G L 2 (0, 1) the opera- 
tors Td(ct) and T^(a) are self-adjoint, bounded below, and have simple discrete spectra. 
We denote by Ai < A2 < • • • the eigenvalues of the operator Td and by p\ < p 2 < • • • 
those of Tn and recall that these eigenvalues interlace, i.e., p n < \ n < p n +i for all 
nGN, and satisfy the relations 

(2.2) y/X^ = nn + p 2 n, y/Ju = 7r(n - \) + p 2n -i 

with some £ 2 -sequence (p n ) [52|l53]. If q belongs to W 2 (0, 1) for some s G [0, 1], then 
by [29] there exists a unique function a* G W 2 (0, 1) such that p n = s n (a*), where 

s n(f) '■= / f(x) sinnnx dx 
Jo 

is the n-th sine Fourier coefficient of /. 

In the paper [51], the mapping ^ sin : / \-t (s n (/))^ =1 defined on the Sobolev 
spaces W 2 (0, 1) was studied in detail for all s > 0. For s G [0, 1] the results of [M] 
can be specified as follows. Denote by W 20 (0, 1) the subspace of W 2 (0, 1) consisting of 
functions that vanish at the endpoints x = and x = 1 and set 

W? j0 (0,l) := [W^ (0,1), 1^(0,1)]^ se (0,1), 

to be the interpolation space [351 Ch. 1.9]. In particular, for s < \ the space W 20 (0, 1) 
coincides with W 2 (0, 1), for s > | it is the proper subspace of the latter consisting of 

1 /2 1 /2 

functions vanishing at the endpoints, and W 20 (0, 1) is a proper subspace of W 2 (0, 1) 
defined by more complicated conditions. We also set 

t 2 := {x = ( lB )^ I ||x|| s 2 := $> 2 '|x n | 2 < 00}; 

this is a Hilbert space under the scalar product (x, y) s := ^2n 2s x n y^, for x := (x n )^ =1 
and y := (y n )'^ =1 . The interpolation theory then shows that for all s G [0, 1] the map- 
ping jF sin is an isomorphism between W 20 (0, 1) and £ 2 ; moreover, under an equivalent 
norm on W 20 (0, 1) this mapping becomes isometric. 

We set Po(x) := 1 — x and P\{x) = x; then for an arbitrary / G W^O, 1) the function 

/„(*) :=f(x)-f(0)P o {x)-f(l)P 1 (x) 

belongs to V^ 2 (0, 1); thus ^(J) G i\ + ls{e , ei}, with 

e, := ^sin(P,) = ((-l^VC™))^ 
and lsS 1 standing for the linear span of a set S. We now set £ 2 := t\ for s G [0, |) and 

^:=^ + ls{e ,e 1 } 
for s G [~, 1]; in this latter case the scalar product in t 2 is introduced via 
(x + a e + aie x ,y + 6 e + Mi) s : = ( x ,y)s + a o&o + ai&i- 
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Then by Lemma 1 of [51] we see that J^sm extends to an isomorphism of VK|(0, 1) 
and t\ for all s G [0,1], which, moreover, is even isometric under an equivalent norm 
on W|(0, 1). Therefore, for q G W£~ (0, 1) the sequence (p n ) defined by (12. 2p belongs 

to 4 

Next we define the norming constants. For a nonzero z G C, we denote by y(-,z) 
a solution of the equation ^(y) = z 2 y satisfying the initial conditions y(0) = and 
yW(0) = 0. Then \An) is an eigenfunction corresponding to the eigenvalue X n of 
the operator T D , and we call the number ot n := -\An)||i2 ^ ne norm ^ n 9 constant 
for this eigenvalue. It is known [28] that if q G W / 2 S_1 (0, 1) for s G [0, 1], then 

Ot n = 2 + /?2n, 

where /?2n is the 2n-th cosine Fourier coefficient C2„(cx) for a unique function a G 
1^1(0, 1) of zero mean that is even with respect to x — |, i.e., cr(l — a;) = cr(x); here 
we set 

c n(f) '■= / f(x) cosnnx dx. 
Jo 

By the arguments similar to the above the mapping J^ cos : / i-)- (c n (/))^ 1 is an iso- 
morphism between the subspace V^ 2 s e ven(0? 1) of even (with respect to x = |) functions 
in W / 2 S (0, 1) of zero mean and the subspace ^ 2even °f sequences in t 2 with vanishing odd 
entries. Clearly, the spaces l\ even and l\ are isomorphic. 

The norming constants a n can be determined from the spectra of the operators Tb 
and Tn as follows. We set S(z) := y(l,z) and C(z) := Due to the integral 

representations (13.11) and (13.21) below, S and C are entire functions of exponential type 1 
with zeros 0, ±a/A^ and 0, iy 7 /^ respectively. The Hadamard canonical products of S 
and C are 

(2.3) s(z) = zJl^- i C(z) = zH 

- LJ - TT n ±J - TT z [n — ±r 

n=l n=l v 

so that 5 1 and C are uniquely determined by their zeros. Finally, we have 

2v A n 



(2.4) a r 



where the dot denotes the derivative in z. 

2.3. Main results. Without loss of generality, we may consider only uniformly pos- 
itive operators, adding the number pi + 1 to q and all the Dirichlet and Dirichlet- 
Neumann eigenvalues as required. Respectively, we introduce the set JV S of data 
((An)neN, (l^n)nen) with the following properties: 

(Nl) pi > 1 and the sequences (A n ) and (p n ) strictly interlace, i.e., p n < X n < p n+ i 
for all neN; 

(N2) the numbers p 2n '■= — nn and p2„_i '■= — 7r(n — |), n G N, form a 

sequence (p„)„ e N in £ 2 . 

In this way every element v := (A, fi) of *yf s is identified with a sequence (p ra ) in l\ or ) 
equivalently, with a unique function / G W% (0, 1) satisfying s n (/) = p n . This induces 
a topology on J\f s \ moreover, ||(p n )||s or ||/||w| define equivalent metrics on jV s . 

According to [28], every element of JY S consists of eigenvalue sequences of the op- 
erators Td(ct) and T^(a) corresponding to some real- valued regularized potential a G 
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W 2 (0, 1) and, conversely, for every real- valued a G W 2 (0, 1) with Xn((x) > I the cor- 
responding spectral data form an element of jY s . When the regularized potential a 
varies over a bounded subset of W 2 (0,1), then the main theorem of the paper [55J 
implies that the corresponding spectral data ((A n ), (/x n )) remain in a bounded subset 
of jV s . Moreover, the Priifer angle technique used in [55J yields then a positive h such 
that all the corresponding spectral data ((A n ), (fi n )) are h-separated, i.e., such that the 
inequalities y/p n +i — > h and \/X^ — sfjl^ > h hold for every n G N. Summarizing, 
we conclude that the uniform stability of the inverse spectral problem we would like 
to establish is only possible on the convex closed sets jV s (h,r) of the spectral data 
consisting of all elements of Jf s that are /i-separated and satisfy ||(p n )||s < r - 
In these notations, one of the main results of the paper reads as follows. 

Theorem 2.1. For every s G [0, 1], h G (0,7r/2) ; and r > 0, the mapping 

^V s (h,r) Bis^ae W^(0,1) 

is analytic and Lipschitz continuous. 

Lipschitz continuity means that there exists a number L = L(s,h,r) such that for 
any two elements Vi and v 2 of ^/V s (h, r) the regularized potentials o~\ and a 2 in W 2 (0, 1) 
solving the inverse spectral problems for the data v\ and v 2 satisfy 

iFl — 0211^1(0,1) < L\v x — VzWjV*- 

See [9] for definitions and properties of analytic mappings between Banach spaces. 

In fact, we prove first analyticity and local Lipschitz continuity in the inverse spectral 
problem of reconstructing a from the spectrum (A n ) of Tj)(a) and the norming constants 
(a n ) (see Theorem 12.21 below), and then derive Theorem 12.11 by showing that the 
norming constants depend analytically and locally Lipschitz continuously on the two 
spectra. 

More exactly, we denote by J?f s the family of strictly increasing sequences A := (A n ) 
for which the sequence (p n ) with p2 n -i = and p 2n '■= \fKi — ^n forms an element of i\ 
and introduce the topology on J?f s by identifying such A with (p n ) G l\. It follows from 
the results of Subsection 12.21 that the subspace ^2,even °f ^2 consisting of elements with 
vanishing odd entries and the subspace W 7 ^ odd(0' 1) °f functions in W^iO, 1) which are 
odd with respect to x — | are isomorphic under the sine Fourier transform J^ S in- For 
h G (0, 7r) and r > 0, we denote by ^f s (h, r) the closed convex subset of Jzf s consisting 
of sequences (A n ) with Ai > 1, \J A n+ i — \f\n > h, and such that ||(p n )||s < r - 

Next, we write £/ s for the set of sequences (a n ) ne N of positive numbers for which the 
sequence (/3 n ) ng psj, with fan-i — and (3 2n '■= cn n — 2, belongs to £ 2 even . This induces the 
topology of £ 2 on £/ s ; by the results of Subsection 12.21 the space ^| even (and thus &/ s ) 

is homeomorphic to the subspace W 2 s even (0, 1) of functions in W 2 (0, 1) of zero mean 
which are even with respect to x = |. We further consider closed subsets £/ s (h,r) of 
stf s consisting of all (a n ) satisfying the inequality a n > h for all n G N and such that 
\\Wn)\U<r. 

It follows from [2S] that to every (A, a) G ££ s x stf s there corresponds a unique 
regularized real-valued potential a G W 2 (0, 1) such that A and ex are the sequences of 
eigenvalues and the corresponding norming constants of the operator Td(ct). The more 
elaborate properties of this mapping are given by the following theorem. 
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Theorem 2.2. For every s G [0, 1], h G (0,7r), h' G (0,2), and any positive r and r', 
the inverse spectral mapping 

^ s {h, r) x £? s (ti, r') 3 (A, a) m- a G W 2 a (0, 1) 

zs analytic and Lipschitz continuous. 

We conclude this section by observing that analyticity and other properties of the 
direct spectral mapping 

W°(0, 1)9(74 (A, ct)e,yV s x 
at least for the classical case s = 1, are well known and studied in detail in, e.g., [4"o] . 

3. Solution of the inverse spectral problem via the GLM equation 

In this section, we recall shortly the method of reconstructing the regularized po- 
tential a based on the Gelfand-Levitan-Marchenko (GLM) equation [25]; see also [51] 
for an alternative approach. Using this method, we shall later show analyticity and 
Lipschitz continuity of the inverse spectral mapping. 

We recall that y(-, z) stands for the solution of the equation l a (y) = z 2 y satisfying the 
initial conditions y(0, z) = and y^(0, z) = z. This function has the representation [27] 

PX 

(3.1) y(x, z) = sin zx + / k(x, t) sin z(l — 2t) dt, 

Jo 

where the kernel k is lower-triangular (i.e., k(x, t) = a.e. in the domain < x < t < 1) 
and has the property that, for every fixed x G [0,1], the functions k(-,x) and k(x, •) 
are in L 2 (0, 1) and depend continuously therein on x G [0, 1]. Also, 



y^\x, z) = z cos zx + z / k 

Jo 



(3.2) y^\x, z) = z cos zx + z / ki(x, t) cos z(l — 2t) dt 



where a kernel k\ has similar properties. Denoting by K the integral operator with ker- 
nel k and by / the identity operator in L 2 (0, 1), we see that I + K is the transformation 
operator mapping solutions of the unperturbed (a = 0) differential equation lo(y) = z 2 y 
to those of l a (y) = z 2 y. 

The GLM equation relates the spectral data for the operator Tb = 7b (er) (i.e., its 
eigenvalues and norming constants) with the transformation operator I + K. To derive 
it, we start with the resolution of identity for Tb, 

N 

I = s-lim V] a n { ■ , y n )y n , 



JV->oo ■ 

71=1 



where s-lim stands for the limit in the strong operator topology, (■,•) is the scalar 
product in L 2 (0, 1), and y n := y(-, y/X^). Recalling that y n = (I + K)s n with s n (x) = 
sin y/X^x, we get 



N 

(3.3) I = (I + K) [s-lim V a „( • , s n )s n ] (I + K*) 



n=l 



The operator in the square brackets has the form I + F, where F = F(X, a.) is an 
integral operator of Hilbert-Schmidt class ©2 with kernel 

(3.4) f{x,t):=l(<t>{x + t)-<j>(\x-t\)), 
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where 

(3.5) <p(x) := (2 cos nnx — cos \[\~ n x) 

ngN 

is a function in L 2 (0, 2). Applying (/ + K*)' 1 to both sides of (13.31) and rewriting the 
resulting relation in terms of the kernels k and /, we get the GLM equation 

/•X 

(3.6) k{x,t)+f(x,t) + h(x,Z)f(Z,t)dt = 0, x>t. 

Jo 

The algorithm of reconstructing q from the spectral data ((A n ), (/x n )) proceeds now 
as follows. We first calculate the numbers a n via (12.41) . then construct the function 
of (13.51) . form the kernel / of (13 .4p . solve the GLM equation (13.61) for k, and set 

rX 

(3.7) a{x) : = -<j>{2x) - 2 / k(x, £)/(£, x) d£. 

Jo 

Then a G L 2 (0,1) is the unique regularized potential for which the operators Tb(cr) 
and Tn(ct) have eigenvalues A n and /i„, nGN, respectively, see [25j[26]. Moreover, if k 
and 4> are smooth, then the GLM equation implies that a(x) = 2k(x,x) — 0(0), thus 
yielding the classical relation 

q{x) = 2-^-k(x, x) 

ax 

for the potential q. It was proved in [28J that if the sequences (A n ) and (fj, n ) are such 
that the corresponding sequence (p n ) is in l s 2l with s G (0, 1], then a of (13.71) belongs 
to W^O, 1), i.e., the reconstructed potential q is in W / 2 S_1 (0, 1). 

4. Reconstruction from a spectrum and norming constants 

In this section, we prove Theorem 12.21 on analyticity and Lipschitz continuity in 
the inverse spectral problem of reconstructing the regularized potentials of Sturm- 
Liouville differential expressions from their Dirichlet spectra and the corresponding 
norming constants. 

We shall study the correspondence between the data (A, a) e ^f s (h,r) x £/ s (h',r') 
and the regularized potentials a G W|(0, 1) of the Sturm-Liouville operator T D (a) 
through the chain of mappings 

(A, cx)\-¥(j)\-¥F\-¥K\-^a, 

in which G ^(0, 2) is the function of (I3.5p . F is the operator with kernel / of (13.41) . 
K G ©2 is the integral operator with kernel k that solves the GLM equation (13. 6p . and, 
finally, a is given by (13.71) . Throughout this section, we fix h G (0,7r), h' G (0,2), and 
positive r and r', and shall always denote by X n the elements of a sequence A G ^f s (h, r) 
and by a n those of an ck G £/ s (h', r'). 

Lemma 4.1. The mapping 

^ s {h, r) x #f s (ti, r') 3 (A, a) m- G W°{0, 2) 

is analytic and Lipschitz continuous. 

Proof. We have 4>{2x) = <£\(x) + ip\, a (x), where 



fx(x) := 2 y^[cos(27rna:) — cos(2a/A^x)] 



n=l 
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and 

oo 

■= - ^/3„cos(2a/A^x), 

n=l 

with (3 n := a n — 2. We recall that a/A„ = irn + S2 n (f) for a unique / G W| odd (0, 1) 

and /3 n = C2 n {g) for a unique g G W-f even (0, 1) an d that the mapping (A, a) >->■ (/, g) is 

isomorphic from «5f s x ,g/ s into W| odd (0, 1) x W2 even (0, 1). Therefore analyticity and 
Lipschitz continuity of the mapping under consideration follows from Corollary IC. 31 □ 

It is advantageous to regard the GLM equation (13. 6p as the relation between the 
integral operators K and F generated by the kernels k and /. We shall need several 
related notions, which we now recall. 

The ideal ©2 of Hilbert-Schmidt operators consists of integral operators with square 
summable kernels, and the scalar product (X,Y)2 '■— tr(XY*) introduces a Hilbert 
space structure on ©2. For an integral operator T G ©2 with kernel t we find that 
(T, T) 2 = Jq Jq 1 \t(x, y)\ 2 dx dy; thus the estimate 



(4.1) 



\f{x,y)\ 2 dx<\ I \<f>( x + y)\ 2 dx+± / \cf>(\x-y\)\ 2 dx 



2 

11 JO 

f2 



< / 10(0 1 2 ^ 
Jo 



shows that F G ©2 and ||-F||g 2 = (F,F) 2 < ||0Hi 2 (o2)- Moreover, the mapping (f> 1— > F 
is linear (and thus analytic and Lipschitz continuous) from L 2 (0,2) into ©2- 

Further, we denote by ©^ the subspace of © 2 consisting of all Hilbert-Schmidt 
operators with lower-triangular kernels. In other words, T G ©2 belongs to ©^ if the 
kernel t of T satisfies t(x, y) — a.e. outside the domain Q + := {(2, y) | < y < x < 1}. 
For an arbitrary T G ©2 with kernel t the cut-off t + of t given by 



t + (x,y) 



t(x, y) for x > y, 
for x < y, 



generates an operator T + G ©J, and the corresponding mapping V + : T i-> T + turns 
out to be an orthoprojector in © 2 onto ©J, i.e. (V + ) 2 = V + and (V + X, Y) 2 = 
(X, V + Y) 2 for any X, Y G © 2 ; see details in [16, Ch. 1. 10]. 
With these notations, the GLM equation (13. 6p can be recast as 

(4.2) K + V + F + V + (KF) = 
or 

(1 + V+)K = -V + F, 

where V\ is the linear operator in © 2 defined by V^Y = V + (YX) and I is the identity 
operator in © 2 . Therefore solvability of the GLM equation is strongly connected with 
the properties of the operator V%. 

Lemma 4.2. For every X G ^(L 2 (0, 1)) ; the operatorV\ is bounded in © 2 . Moreover, 
for every F from the set 

(4.3) F := {F(X, a) | (A, a) G ^ s {h, r) x ^ s (h', r')} 

the operator X + Vp is invertible in ^(©2) and the inverse {X + Vp)~ l depends ana- 
lytically and Lipschitz continuously in £${f5\) on F G F in the topology of G 2 - 
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Proof. Boundedness of V x in ©2 is a straightforward consequence of the inequality 

||^X*1I©2 - ll^^ll©2 - ||-X"||^(i 2 (0,l))ll^ r ||s25 

see [151 Ch. 3]. Assume next that / + X > el in L 2 (0, 1); then, for any Y G & 2 , 

((X + P+)Y, Y) 2 = (Y, Y} 2 + (YX, Y) 2 = tr(K(J + X)Y*) . 
We see that Y(I + X)Y* > eYY* and by monotonicity of the trace we get 

((X + V x )Y,Y) 2 >e(Y,Y) 2 , 

i.e., X + V x > eX in ©£. 

Now, if F = F(A, a) is constructed as explained in Section [3] from A = (Xk)ken G 
Jz? s (h,r) and a = (afc)fceN 6 £/ s (h',r'), then 

V 

I + -F = s-lim > ajfc( • , s fc )s fc 

fc=l 

with Sjt(x) := sin \f\~kX. By definition, a& > /i' for all fcGN; moreover, by Theorem lA.il 
the sequence (sk)km forms a Riesz basis of the space L 2 (0, 1) and its lower bound m > 
can be chosen the same for all A G JZ' s (h,r). Therefore 

oo 

((I + F)y,y) = ^a k \(y,s h )\ 2 > h'm\\y\\ 2 

k=l 

for every y G £2(0, 1), so that I + F > h'ml. 

By the above, X + Vp > h'mX in (5^; thus X+Vp is boundedly invertible in ^(©2 ) 
and 

||(X + P+)- 1 || < (h'm)- 1 . 
Since V x depends linearly on X, it follows that the mapping 

F^ (Z + P+)- 1 

from ©2 into 3S{& 2 ) is analytic and Lipschitz continuous on the set F. The proof is 
complete. □ 

Corollary 4.3. For every F G J 7 ', t/ie GLM equation ( 14. 2 p /ias a unique solution 

K := -(X + V^Y X V+F G 6+; 

moreover, the mapping F \-t K from F C & 2 to & 2 is analytic and Lipschitz continu- 
ous. 

In view of the above results and formula ( 13. 7ft . the regularized potential a is deter- 
mined uniquely by the function of (13. 5p . To complete the proof of Theorem 12.21 we 
shall show that the induced mapping <fi 1— )■ a is analytic and locally Lipschitz continuous 
from the space W 2 S (0,2) into W 2 S {0, 1) for every s G [0, 1]. We shall establish this for 
s = and s = 1, and then interpolate to cover all the intermediate values s G (0, 1). 

4.1. The case s = 0. This case is the easiest to treat, although it corresponds to 
the set of Sturm-Liouville operators with the most singular potentials — namely, with 
distributional potentials in W / 2 " 1 (0, 1). 

Lemma 4.4. The function a of ( 13. 7p depends analytically and Lipschitz continuously 
in X 2 (0, 1) on the function (p G L 2 (0, 2) of ( 13. 5 p that is in the range of the mapping of 
Lemma \4-l\ for s = . 
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Proof. By definition, a(x) = — 0(2x) —2 J* k(x, t)f(t, x) dt, where k is the kernel of the 
solution K of the GLM equation ( 13. 6 j) and f(t, x) = \[4>{t + x) — (f>(\x — t|)]. Thus the 
integral above depends linearly on k and 0; moreover, in view of (14.1 j) we get 

•1 /-a 2 

fc(x, £)/(£, x) (it 



2 J + / I j/j. „\|2 



fix 



< dx \k(x,t)\ 2 dt / |/(i,x)| 2 dt 
Jo </o Jo 

< 110111(0,2) / / \k(x,t)\ 2 dxdt= ||0||! 2(O)2) ||K||| 2 . 

Since if G depends analytically and Lipschitz continuously on F G F and F G ©2 
depends linearly and continuously on = 0(A, a), the result follows. □ 

This completes the proof of Theorem 12.21 for the case s = 0. 

4.2. The case s = 1. The function of (13. 5p belongs in this case to IV^O, 2), and the 
solution k to the GLM equation (I3.6P must also possess some extra smoothness. We 
recall that functions in W^O, 1) are continuous and that there is C > such that 

max \g(x)\ < C||#||wi (0 1} 

»e[o,i] 2 

for every g G W^O, 1). As above, we denote by Q + the set {(x,t) | < t < x < 1}. 

Lemma 4.5. Let s = 1 and and k be defined as above. Then the distributional 
derivative d x k of k belongs to L 2 (fl + ) and the induced mapping (A, a) h-> d x k is analytic 
and Lipschitz continuous from ^' 1 (h,r) x ^/ 1 (h',r') into L 2 (fi + ). 

Proof. Let (A, a) be the element of J§f r) x &/ l {h' , r') that generates the function 
G W 2 X (0,2). We set 



[X) := 

fc=i 



(2 cos 7r/;;x — cos a/A^x) , 



which corresponds to taking = ctk,o '■= 2 and A& = A^o := 7r 2 /c 2 for k > n. Choose 
n so large that \y/X^ — irn\ < n — h if n > n ; then for such n the sequences 

(A*)Li U (A k,0)k>n 

and 

(0ik)k=l U (ttfc,o)fc>n 

belong to J£ 1 (h, r) and ! g/ 1 (/i',r / ) respectively and n — >■ in W^O, 1) by Lemma [4. II 
We form the kernel /„ taking n instead of in (j3.4P ; then the GLM equation with / 
replaced by /„ has a unique solution k n G L 2 (fi + ). Denoting by F n and if n the integral 
operators with kernels f n and k n , we see that / n — > f in AmO, 1) x (0, 1)) means that 
F n — >■ F in & 2 ; hence K n -f K m. & 2 by Corollary I4T3| i.e., fe„ ->• k in L 2 (fi + ). 

Since the integral operator F n is of finite rank, the solution k n can be written in 
an explicit form and is easily seen to be smooth in the domain Q + , cf. [TU Sect. 12] 
and [HI Ch. IV. 3]. We set l n := d x k n ; then l n satisfies in Q + the equation 

l n (x,t)+ / l n (x,£)f n (£,t)d£= -f n (x,t) - k n (x,x)f n (x,t), 

Jo 
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where f n := d x f n . The convergence n — > in W^O, 2) implies that f n — > f in C(n + ) 
and f n (x,t) -> := |[0'(a; + t) - ^'(a; - t)] in L 2 (fJ + ); also, 

a n (x) := 2/c n (x, x) - n (O) = -0„(2x) - 2 / fc n (x, t)f n (t, x) dt 

Jo 

converge in L 2 (0, 1) to cr(x) by Lemma [4.41 It follows that the kernels 

g n (x, t) := f n (x, t) + k n (x, x)f n (x, t) 
converge in L 2 (f2 + ), as n — > oo, to 

g{x,t) := f(x,t) + \[a{x) + <f>(0)]f(x,t), 

with / := d x f. Denoting by L n , G n , and G the integral operators in <5^~ with kernels 
l n , g n , and g respectively, we conclude by Lemma I4T21 that 

L n = -(X + P+ Y l G n -)• -(X + P+)- 1 G =: L 

as n — > oo in the topology of Therefore l n converge in L2(fl + ) to the kernel / 
of the operator L. We conclude that, in the sense of generalized functions, / = d x k 
and d x k belongs to L 2 (fi + ) as claimed. It is easily seen that the mapping (A, a) y G 
from = Sf 1 (/i, r) x &/ 1 (h' , r') to & 2 is analytic and Lipschitz continuous. We finally apply 
Lemma 14.21 to conclude that L G © 2 depends in the same manner on F G T and 
G G &2', here T is the set of (14.31) corresponding to s = 1. □ 

To complete the proof of Theorem 12.21 for s = 1, it suffices to show that the function 
o"i(x) := JJf k(x, t)f(t, x) dt depends analytically and Lipschitz continuously in W^i^O, 1) 
on (A, ot) G =$f 1 (/i, r) x £^ 1 {h\r l ). That o"i depends in this manner in L 2 (0, 1) was 
established in Subsection 14. 11 Also, 

/•X PX 

= WW + 0(O)][0(2x) -0(0)] + / l(x,t)f(t,x)dt+ / k(x,t)f(t,x)dt, 

Jo Jo 

where, as in the proof of Lemma 1431 l(x,t) := d x k(x,t) and f(t,x) := d x f(t,x). 
Clearly, the first summand above belongs to L 2 (0, 1) and depends therein analytically 
and Lipschitz continuously on (A, a). Also, / and / depend in the same manner 
in L 2 (f2 + ) on the spectral data (the former by Lemma [4.5[ the latter by linearity and 
direct estimates (14. ip ). Thus both integrals give functions in L 2 (0, 1) with required 
dependence on (A, a.) (see the proof of Lemma T4. 41) . which establishes Theorem 12.21 for 
s = 1. 

4.3. The case s G (0, 1). The general case will be covered by the nonlinear interpola- 
tion theorem due to Tartar |58j, which implies the following result. 

Proposition 4.6. Assume that (X , Xi) and (Y^^i) ore pairs of Banach spaces with 
continuous embeddings X± > X and Yi •—^■Yq. Let also $ : X — >• Y be a nonlinear 
mapping that is Lipschitz continuous on the balls B Xo {r) := {x G X \ \\x\\ Xo < r} for 
every r > 0. Assume further that QXi C Yi and that $ is Lipschitz continuous on 
every ball Bx 1 {r) of X\ as a mapping from X\ into Y±. Construct the interpolation 
spaces X s := [Xi,A"o] s and Y s := \Yi,Yq] 8 , s G (0,1), by the complex interpolation 
method; then $ acts boundedly from X s to Y s for every s G (0, 1) and, moreover, its 
restriction to the ball Bx s {r) of X s is Lipschitz continuous for every r > 0. In other 
words, for every r > there is C = C (r, s) such that 

mx 1 )-Q{x 2 )\\ Ym <C\\x 1 -x 2 \\ x . 
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whenever x\ and X2 belong to Bx B {r). 

By definition, the spaces l\, t\ and W|(0, 1), as well as their "even" and "odd" 
subspaces, form the interpolation space scales. Therefore the above proposition, in view 
of the results of Subsections 14.11 and I4.2[ implies that for every h G (0,7r), hi G (0,2), 
positive r and r', and s G (0, 1) the mapping (A, a) i— > a is Lipschitz continuous 
from J5? s (h,r) x £/ s (h',r') to W|(0, -0- Analyticity for all s G [0, 1] again follows from 
that for s = and s = 1 and interpolation theorem for linear operators [6,35j applied 
to the Frechet derivative of this mapping. 

Remark 4.7. For s < ~, the arguments of the paper [2Sj combined with the results 
of the previous subsection provide a direct proof of analytic and Lipschitz continuous 
dependence of a G W^O, 1) on the spectral data (A, a) G J£ s (h,r) x g/ s (h',r'). 

Remark 4.8. It should be clear how one can iterate the considerations of Subsection 14.21 
to all natural values of s and then interpolate as in Subsection 14.31 to get all positive s, 
cf. [381 Sect. 3.4] and [57]. 

5. Reconstruction from two spectra 

In this section, we complete the proof of Theorem 12.11 by establishing uniform con- 
tinuity of the norming constants on the two spectra. More exactly, given an ele- 
ment ((A ra ), (yU n )) G ,jV s [h, r), we define numbers a n by the relation 

(5.1) a r , " X ; "' 



where the entire functions 5* and C are given by the infinite products (I2.3p . and prove 
the following theorem. 

Theorem 5.1. For every s G [0, 1], h G (0,7r/2) and r > the mapping 

^ s {h,r) 3 (A,/x) a G srf s 

is analytic and Lipschitz continuous; moreover, there are h! > and r' > such that 
the range of this mapping belongs to £/ s (h',r'). 

It suffices to show that there are h" > and r" > such that the numbers (3 n defined 
via the relation 

2 

— = 1 + 0n 

satisfy the inequality 1 + j3 n > h" and form an £|-sequence (3 := ($ n ) depending 
analytically and Lipschitz continuously on (A, fi) G ^Y s (h,r) and having norm not 
greater than r". Indeed, then a n > h! := 2/(1 + ?-") and /?2n := «„- 2 = — 2/3 n /(l+/3 n ). 
We observe that, for a bounded sequence (d n ), the mapping (x n ) !->■ (d n x n ) is a bounded 
linear operator in t\ of norm d := sup n \d n \. Since |— 2/(1 + /5 n ) < 2//i", we conclude 
that the sequence (/?2n) belongs to £| and has the norm at most 2r"/h"; the sequence 
(j3 n ) with /?2n-i — belongs then to ^| even and is of norm at most r' := 2 s l 2+1 r" /h" . 

To justify analyticity and Lipschitz continuity of o: in &/ s , we exploit the fact that 
t\ is a Banach algebra under the point-wise multiplication (x n ) ■ (y n ) = (x n y n ). We 
denote by A the unital extension of l\] elements of A have the form al + x, where 
a G C, 1 is the sequence with all its elements equal to 1, and x = (x n ) G £|, and the 
norm in A is given by ||al + x]]^ := |a| + ||x|| s . An element al + x is invertible in A 
if and only if a ^ and a + x n ^ for all n G N; in this case the inverse is equal to 
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-1 + y with y = (y n ) and y n := —x n /[a(a + x n ]. Since the sequence (l/[a(a + x n )]) is 
bounded, the above reasoning shows that indeed y belongs to l\. 
We now see that 

ia = (l + /3)- 1 ; 

since taking an inverse element is an analytic mapping in a unital Banach algebra, a 
depends analytically on (3. Lipschitz continuity of ct follows from the fact that, for the 
set of f3 considered, 1 + (3 have uniformly bounded inverses in A (of norm not greater 
than 1 + r"/h"). 

To the rest of this section, we shall use the notations u 2n '■= \fKi and U2 n -i '■= y/JJm 
and p n := u n — nn/2. In view of ( 15. ip we have 

1 +/3 n = S(u 2n ) =: (1 + On) (1 + b n ), 

where we set 

(5.2) a n := (-l) n S(co 2n ) - 1, b n := (-1)-^^ - 1. 

Thus we need to prove that both S{uj2 n ) and C{ui2 n ) /u 2n are uniformly bounded away 
from zero, that the sequences (a n ) and (b n ) are the sequences of the even cosine Fourier 
coefficients of some functions hi and h,2 from W 2 eve n^i 1) °f zero m ean, and that the 
mappings 

(5.3) J^ s (h, r) 3 (A, n)^h x e Wl even (0, 1) 
and 

(5.4) jV s % r) 3 (A, M ) m- h 2 E ^, even (0, 1) 

are analytic and Lipschitz continuous. We do this in the two subsections that follow. 
5.1. Analyticity and continuity. 

Lemma 5.2. The mappings of ( 15. 3 j) and (15. 4p are analytic and Lipschitz continuous. 

Proof. We observe that in view of (13. ip and (13. 2p the functions S(z) and C(z)/z have 
the representation 

cos z+ / g{t) cos z{\ — It) dt. 



with g(t) = (l-2t)k{l,t) E L 2 (0, 1) for S(z) and g(t) = h(l,t) E L 2 (0, 1) for C{z)/z. 



Therefore both a n and b n can be written as 



[(-l) n cosw 2n - 1] + (-1)™ / g(t) cos u 2 n(l-2t)dt=:d n + e n , 

Jo 

with respective g E L 2 (0, 1). Clearly, we may (and shall) take the even part g eve n of g 
instead of g in the above integral. 

Recalling that u 2n = 7m + s 2n (f) for a (unique) function / E W2 odd (0, 1) and ob- 
serving that s 2n (f) = with f{n) being the n-th Fourier coefficient of a function / 
(see Appendix |B]), we can write the d n as 

°° f(n) 2k 

d n = cos s 2n (f) - 1 = cosh f(n) - 1 = . 

fc=i 
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It follows that the d n is the n-th Fourier coefficient of the function / given by 

j<2k> 



; (2*)! ' 

where f <k> is the /c-fold convolution of / with itself. The function / has zero mean 
and is even with respect to x — 5 and therefore d n = C2n(f)- By the results of 
Appendix IB1 the mapping / i->- / is analytic from W 2 s odd (0, 1) to W2 evea (0,l) and is 
Lipschitz continuous on bounded subsets of /. 

Recalling Lemma IC.4I and the remark with formula (IC.2I) following it, we see that 
the e n give the 2n-th cosine Fourier coefficient of the function h := §\fr(i/, g) in 
W| even (0, 1) °f zero mean, which depends analytically and boundedly Lipschitz contin- 
uously on / and g. 

It remains to apply Lemma [5731 below to conclude that the even parts of the functions 
(1 — 2t)k(l,t) and ki(l,t), which we have taken as the function g above, depend 
analytically and Lipschitz continuously in W 2 S {0, 1) on the spectral data in ^Y s (h,r). 
This completes the proof of the lemma. □ 

Lemma 5.3. For every s G [0, 1], the mappings 

Jf s {h, r) 3 (A, ft) ^ fc od d(l, •) e W 2 S (0, 1) 

and 

jK s (h,r) 3 (A, fi) ^ k hevcn (l r ) G W 2 S (0,1) 

are analytic and Lipschitz continuous. Here fc dd(l, •) and fci je ven(lj •) ar ^ respectively 
odd and even parts (with respect to x = \) of the functions k(l, •) and k±(l, •). 

Proof. Since both mappings are treated similarly, we only consider in detail the first 
one. By definition of the function y(-,z) we have y(l,u) 2 n) — 0, and thus the numbers 
^2n = vrn + s 2n (/) along with and — a> 2n are zeros of the entire function 

y(l,z) = smz+ [ k(l,t)sxnz(l-2t)dt. 
Jo 

This function does not depend on the even part k even (l, •) of the function k(l, •); in 
fact, it can be written in the form 

(5.5) sinz+ I g(t)e lz{l - 2t) dt 

Jo 

with g{t) := —\k Q ^{l,t). Since the set of functions {sina; 2n (l — 2i)} neN is complete 
in ^2,odd(0, 1) (cf. the results of Appendix IA1 . one can show that fc dd(l, •) is uniquely 
determined by the zeros u 2n = vrn + s 2n (/) = irn + if in). Some important properties 
of the induced mapping / h-> g can be derived from the paper [3D] . 

Indeed, the results of [30] imply that for every / G W^O, 1) there exists a unique 
function g G Wf(0, 1) such that all zeros (counting multiplicities) of the entire func- 
tion (15.51) are given by the numbers rcn + f(n), n G Z. Such pairs of / and g in fact 
satisfy the relation 

(5.6) H(f, g) := s(f) + g + £ W*/^ = 0; 

fc=i 
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here 

00 t 1 \k f<2k+l> 

s (f) ■= y izlj I 

S[Jh 2- (2k + l)\ ' 

/ <fc> is the fc-fold convolution of / with itself, and M is the operator of multiplication 
by i(l - 2x). The function H is analytic from W|(0, 1) x t^ 2 s (0, 1) into W^O, 1), and 
its partial derivatives dfH(f,g) and d g H(f,g) are given by 



(5.7) fylftf, g){h t ) = (c(f) + £ 1 9 >J ) * hi 



k=l 

00 

(5.8) d g H(f,g)(h 2 ) = h 2 + J2 

k=l 



(M k g)*f <k - 1 
{M k h 2 ) * f <k> 



with 



if) ■■= E 



k £<2k> 



i)7 



,..„ < 2 *> ! 

We assume now that / £ W 2 s odd (0, 1) * s suc h that the corresponding sequence A = 
(A n )neN with A n := u 2n = + s 2n (f)) 2 belongs to J£ s (h, r). Set 5a to be the function 
of ( 12. 3ft ; then S\ can also be represented as (15. 5p . Direct calculations show that the 
n-th Fourier coefficient of the function of ( 15 . 71) is equal to 

{-l) n k{n) \cosu 2n + / i(l - 2t)g(t)e^ 1 ^ dt] = (-^(n^KO- 







By Lemma 15.41 below there are numbers K\ and K 2 such that 

^i<|5 A (u; 2n )|<^2 

for all A £ =£f s (/i, r) and all neN. The results of Appendix |B] imply that the operator 
dfH(f,g) is bounded in every space W 2 (0, 1) and its norm is at most K 2 . 
Similarly, the n-th Fourier coefficient of the function of ( 15. 8 p is equal to 

(-l) n I h 2 (t)e iul2n(1 - 2t) dt. 



By Theorem IA.1I there exist positive M and m such that, for all A £ J2? s (/i, r), the 
sequences {<d lU2n ^ l ~ 2x ') n & form Riesz bases of £2(0, 1) of upper bound M and lower 
bound m, see Appendix lAl Therefore the operator H g := d g H(f,g), 



H g : h 2 ^ ^(-l) n (h 2 ,^ i2x - 1) )e* nrdx 



is bounded and boundedly invertible in L 2 (0,1), with \\H g \\ < M 1 ! 2 and 1 || < 
mT 1 ! 2 . If h 2 £ W^O, 1), then we integrate by parts to get 

Cn : = (h 2 ,e iu ^ 2x -V) = J—[h 2 (0)^ - h 2 {l)e-^} + - J_(^,e^ 2a, - 1) ). 

2iu 2n 2iu 2n 

It is clear that the sequence (c n ) ng ^ forms an element of £l(Z), see Section IB~2l Thus 
the operator H g acts boundedly and boundedly invertible in W 2 (0, 1), and it remains 
to use the interpolation theorem to derive the same properties of H g in W 2 (0, 1) for all 
s£[0,l]. 

We now use the implicit mapping theorem to conclude that the mapping / 1— > g is 
analytic in W|(0, !)• Recalling the isomorphism of the space =Sf s of sequences A = (A n ) 
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of the Dirichlet eigenvalues of the Sturm-Liouville operators T(a) with a G W|(0, 1) 
and the subspace W^|odd(0' -0 explained in Section [21 we conclude that the mapping 

Sf'(h, r)9A4 fc odd (l, ■) G W? )Odd (0, 1) 

is analytic. The uniform bounds on the inverses of the partial derivatives dfH(f, g) 
and d g H(f,g) established above imply that this mapping is Lipschitz continuous, and 
the proof for the first mapping is complete. 

The second mapping of the lemma is treated analogously using the relations 

coscj2n-i + / ki(l, t) coso;2 n -i(l — 2£) dt = 

JO 

and the Riesz basis properties of the system (cosco^n-i^neN, cf. Remark I A. 2 1 □ 

5.2. Uniform positivity of a n . Since jY s {h, r) C ^V{h, r) if s > 0, it only suffices to 
consider the case s = 0. In view of formula (15.11) . uniform positivity of a n immediately 
follows from the lemma below. 

Lemma 5.4. For every h G (0,7i"/2) and r > we have 

sup sup log IS^c^n)! < oo, sup sup log < 00) 

(A,fi) ™eN (A, M ) neN ^2n 

where S and C are constructed via (I2.3P from sequences A := (oJ% k ) and fi := (w. 
and the suprema are taken over (A, n) G ^V(h,r). 

Proof. By (12.31) . we have 



2k-l)> 



Dividing both sides by 



dsinz 



cos 7m 



dz 

we conclude thatQ 



2 n 



TT 2 k 2 — ir 2 n 2 



7T 2 k 2 
fcSN, k^n 



I^VV^nJI - -2^2 11 I2L2 11 



7r 2 n 2 7r 2 n 2 — 7r 2 A; 2 7r(&; — n) ' 

fc^n fceZ, k^n v 7 

where we set u_k := — for fceN and cuo := 0. Set also (recall that p^ := Uk — irk) 

_ ^2k — ^2n 1 P2fc — P2n 

a fc,n ; — — n \ — 1 



7r(k — n) 7r(k — n) 

ii k ^ n and a nin := 0; then ^(a^n)! = ELez^ + a k,n)- Since the sequence (co n ) is 
/i-separated for every (A, /x) G ^/K(h,r), we have 1 + <2fc in > 2/i/7r for all integer 
and n. Therefore, with 



log(l + x) — X 

K := max 

x>—l+2h/n 



X 2 



< OO, 



we get the estimate 

(5-9) |logJJ(l + a fe , n ) < \j2 a k,n +Kj2 a lm 



fcez fegz fcez 



1 In what follows, all summations and multiplications over the index set Z will be taken in the 
principal value sense and the symbol V.p. will be omitted. 
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provided the two series converge. 
Clearly, 



E- 



k^n 



k — rt 



and thus 



P2fc 



fcez 



7r f— k — n 

kjtn 



< 



V2r 



by the Cauchy-Bunyakovski-Schwarz inequality (recall that J2kezP2k — ^ r ' 2 ^y the 



definition of the set ,jV{h,r) and ^2 k ^ n (k 



n) 



7r /3). Next, the inequality 



4 fc,n 



< 



2p 2 



2A- 



(A; — n) : 



+ 



2p 



2« 



for 7^ n yields 



(A; — n) 2 
1 4vr 2 r 2 



ke 



It follows from (15. 9 j) that 

I log + °*>») ^ + 4AW)/3, 



where the constant only depends on /i. 
Similarly, we find that 



In 



n 



^2fe-l — ^2n 



7T71 



fc=l " v ~ 2/ fcgZ 

and then mimic the above reasoning to establish the other uniform bound. The lemma 
is proved. □ 

As explained at the beginning of this Section, the above statements complete the 
proofs of Theorems 15.11 and 12.11 
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Appendix A. Riesz bases of sines and cosines 

We recall (see, e.g., [151 Ch. 6] and [591 Ch. 4]) that a sequence (f n )nen m a separable 
Hilbert space H is called a Riesz basis of H if it is a homeomorphic image of an 
orthonormal basis of H. Then there are M > (the upper bound) and m > (the 
lower bound) such that, for every / 6 H, we have 



711 



< 



£!(/,/« 



< M 



Riesz bases of L 2 (0, 1) that are composed of exponential functions, or sines, or cosines, 
have been extensively studied in the literature starting from the early 1930-ies, see 
the books by Paley and Wiener [15] and Avdonin and Ivanov [5] for particulars and 
historical comments. For instance, the famous Kadets ^-theorem [3T| implies that for 
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every L < \ there exist positive constants m and M such that as long as a sequence 
(oj n ) n( zz of real numbers satisfies the condition 

(A.l) sup \cj n — nn\ < tyL, 

then the sequence (e luJnX ) n€ z of exponentials forms a Riesz basis of L 2 (— 1, 1) of upper 
bound M and lower bound m. Analogous results for families of sines and cosines were 
established in [18] . 

In this paper, we need generalizations of these results to sequences that may not 
satisfy condition flA.lj) . Recall that Jzf°(/i, r), with h G (0, n) and r > 0, stands for 
the set of all strictly increasing sequences A = (w^) ng N of positive numbers satisfying 
the conditions Ui > h, u) n+ \ — u n > h, n G N, and \ u n ~ nn \ 2 — 7-2 ■ F° r every 
A G Jz?°(/?,,r), we denote by and the sequences of functions (smu} n x) ne w and 
(cosoj n x) ne z+ respectively, with u := 0. We also set w_ n := — u n and denote by S'x 
the sequence of functions (e ltJ "^ 1_2x ^) ne z. The following statement can be derived from 
the results of E 



Theorem A.l. For every h G (0,7r) and r > there exist positive numbers M and 
m such that for every A G J2?°(/i, r) the sequences 5^\, and $\ are Riesz bases of 
£2(0, 1) of upper bound M and lower bound m. 

Remark A. 2. For h G (0,7r) and r > 0, we denote by ^°(h,r) the set of increasing 
sequences /J, := (^n)^Li with the following properties: 
(Ml) Hi > 1 and, for all n G N, y/[i n +i — y/J^n > h; 

(M2) the numbers p n := ^/fI^ — ir(n — |) form a sequence in £2 of norm at most r. 
Then an analogue of the above theorem holds true for the family of sequences ^ = 
(cos v //I^x) ne N, with [i running through the set ^°(h,r); see [2] 



Appendix B. Sobolev spaces W 2 s (0, 1) and some of their properties 

We recall here some facts about the Sobolev spaces W 2 {0, 1) and Fourier coefficients 
of functions from these spaces. For details, we refer the reader to [351 Ch. 1]. 

B.l. The definition. By definition, the space W/^O, 1) coincides with L 2 (0, 1) and 
the norm || • || in W 2 (0, 1) is just the L 2 (0, l)-norm. For a natural I, the Sobolev space 
^(0, 1) consists of all functions / in L 2 (0, 1), whose distributional derivatives for 
k — 1, . . . , I also fall into £2(0, 1). Being endowed with the norm 

(B.l) ll/lk:=(Ell/ (fc) Ho) 1/2 ' 

k=0 

the space W 7 ^^ 1) becomes a Hilbert space. 

The intermediate spaces W 2 (0, 1) for arbitrary positive s can be constructed by 
interpolation [351 Ch. 1.2.1]. We shall need such spaces only for s G [0,2] and thus 
interpolate between W 2 (0, 1) and W-^ ^ 1) to get them, i.e., 

Wf(0, 1) := [Wf (0, 1), W^O, l)]i-t, t G (0, 1). 

The induced norms || • || s (for s = 1 the norm (IB.lj) is equivalent to that defined by 



interpolation) are nondecreasing with s G [0,2], i.e., if s < r and / G W 2 (0,1), then 
\\f\\s < \\f\\r- Since by construction the spaces W2 S (0, 1) form an interpolation scale, the 
general interpolation theorem [221 Theorem 1.5.1] implies the following interpolation 
property for operators in these spaces. 
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Proposition B.l. Assume that an operatorT acts boundedly in W^(0, 1) and W%(0, 1), 
s < r. Then T is a bounded operator in W / 2 * s+( ' 1 ^ (0, 1) for every t G [0, 1]; moreover, 

W 1 \\ts+(l-t)r S || J Hsll-t || r ■ 

Proposition IB. II yields boundedness in every VF|(0, 1), s G [0,2], of the reflection 
operator R given by Rf(x) = f(l — x) and the operator M of multiplication by ix, 
Mf(x):=ixf(x). 

B.2. Fourier transform. For an arbitrary / G ^(0, 1) we denote by / : Z — > C its 
discrete Fourier transform, viz. 

/(n) := I f{x)e" 2mnx dx. 
Jo 

The Fourier transform is a unitary mapping between £2(0, 1) and £ 2 (Z). We set 

f ^|(Z) if s < |, 

§(Z):=^(Z) + ls{e«} if |<«<§, 

|/|(Z) -i-ls{e«,e (2) }, if |<s<2, 

where e^, j = 1,2, denotes the sequence (e^) ne z with := n"- 7 for n 7^ and 
:= 0. The norm || ■ || s in £|(Z) is defined as follows: given an element x := (x n ) 
of £|(Z) and complex numbers a\ and a 2 (with ai = a 2 = if s < | and a 2 = if 

I < s < |), we set 

||x + axe (1) + a 2 e (2) ||2 ■= ^(1 + n 2 ) s |x„| 2 + \a x \ 2 + |a 2 | 2 . 

It is straightforward to verify that the Fourier transform of every function in W$ (0, 1) 
forms an element of ^(Z) and that ||/|| 2 introduces a norm on 1F 2 2 (0, 1) that is equiv- 
alent to that of (IB.ip . Since the spaces ^(Z), s G [0,2], form an interpolation scale, 
the interpolation theorem implies that, for every s G [0,2], the Fourier transform is a 
homeomorphism between the spaces W£(0, 1) and ^f(Z). 

B.3. Convolution. As usual, * denotes the convolution operation on (0, 1), viz. 

(f*g)(x) := f f(x-t)g(t)dt; 
Jo 

here we extend a function / onto the interval (—1,0) by periodicity, i.e., by setting 
f(x) := f(x + 1) for x G (—1,0). It is well known that convolution accumulates 
smoothness; the precise meaning of this statement is as follows. 

Proposition B.2. Assume that s,t G [0,1] and that f G W 2 S (0, 1) and g G W 2 *(0, 1) 
are arbitrary. Then the function h := f * g belongs to H / 2 s+t (0, 1) and, moreover, there 
exists C > independent of f and g such that \\h\\ s+t < C\\f\\ s \\g\\ t . 

Proof of this proposition is based on interpolation between the extreme cases s,t — 
0, 1, which are handled with directly using the representation 

PX rl 

(f*g)(x)= f(x-t)g(t)dt+ f(x-t + l)g(t)dt. 

JO Jx 

We also recall the relation / * g(n) = f{n)g{n). 
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Appendix C. Some auxiliary results 
Lemma C.l. For f and g in £2(0, 1)? set 

Hf,9) ■= V.p.^ff(n)[exp{/(n)ix} - l]e 2 ~ 

then for every s,t G [0, 1] £/ie mapping 

$ : W£(0, 1) x W£(0, 1) ->■ W£ + *(0, 1) 

analytic and Lipschitz continuous on bounded subsets. 

Proof. If / and g are in L 2 (0, 1), then the series for $ converges absolutely, and thus 
$(f,g) is a continuous function. Take now / G VFKO, 1) and g G W^O, 1). Developing 
exp{/(n)ix} into the Taylor series and changing the summation order (which is allowed 
since the resulting double series converges absolutely), we get 

^9) = f:^V.p.J2g(n)f k (n)e^ = £ ^ /*(*), 

k=l ' n&L k=l 

where hf. := g * f <k> , f <1> '■= f, and f <k> for k > 2 is the fc-fold convolution of / 
with itself. By Proposition \B.2\ the functions belong to W 2 +t (0, 1), and 

||^|| s+ *<C fc+1 ||/||^||^||, 

with the constant C of that proposition. Also, the operator M of multiplication by ix is 
continuous in every (0) 1)' r e [0> 2]; denoting by ||M|| r the norm of the operator M 
in W 2 r (0, 1), we conclude that belongs to W 2 +t (0, 1) and 

IW,*)IU < C||^|| t (exp{C||M|| s+t ||/|| s } - 1). 

Moreover, every summand (ix) k hk/k\ depends analytically on / and g; since the series 
converges absolutely in W 2 +t (0, 1), analyticity of $ follows. 
By similar arguments we also find that 

Hh,g) - Hf*g) = E - / 2 <fc> ); 

fc=i 

since 

ll/l <fc> - / 2 </£> lls < - hUWhWs + ||/ 2 ||.)*~\ 

it follows that 

\Mfi,g) - $(f2,g)\\s+t < C^HMIU+tH/i - / 2 |Mb||tex P {C||M|| s+t (||A|| s + ||/ 2 || 5 )}. 
We observe now that the mapping $ is linear in the second argument; therefore, 

- $(f 2 , g 2 )\\ s +t < \Wi,9i) - ®{f 2 , gi)\\ s +t + \Mf2,9i -92)\\.+t 
< CiiWfi - f 2 \\ s + \\ gi - g 2 \\ t }, 

where d < L7 2 ||M|| s+i (l + K) exp{C\\M\\ s+t K} as long as 

\\fi\\s + \\f 2 \\s + \\9i\\t + \\92\\t < K, 
and the desired Lipschitz continuity of $ follows. □ 
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We observe that 

Hf, g)+9 = V.p. 9( n ) exp{/(n)ix}e 2 



and hence the mapping 



(f,g) H- V.p. J2g(n) exp{f(n)ix}e 2 ™ 



is uniformly continuous from W|(0, 1) x W|(0, 1) into W 2 *(0, 1). 

Also, in the definition of $, we can formally take g to be the unity 5 of the convolution 
algebra L2(0, 1), which results in g(n) = 1. Slightly adapting the above proof, we get 
the following result. 

Lemma C.2. Fix s e [0, 1] and, for f e W£(0, 1)> srf 

Kf) := V.p. ^[exp{/(n)ia;} - l]e 27rira . 

TTien t/ie function h(f) belongs to Wf(0, 1) ^ e mapping f i-> zs analytic 

in W£(0, 1) and Lipschitz continuous on bounded subsets. 

We say that a function / is ode? (resp. even) over (0, 1) (with respect to |) if it 
satisfies the relation f(l — x) = —f(x) (resp., the relation f(l—x) = f(x)) a.e. on (0, 1). 
For every integrable function / we define its odd part / Q dd and even part / cvcn by the 
equalities 

_ f( x )-f(l- x ) _ + 

Jodd^J • 2 ' JeveayX) ■ 

Since the reflection operator is continuous in every space W|(0, -Q> s — 0' the mappings 
/ ^ /odd and / i — y f cven are bounded in every W£(0, 1). We denote by W / 2 s odd (0, 1) and 
W| even (0, 1) the subspaces of W|(0, 1) consisting of functions that are respectively odd 
and even over (0, 1). 

Corollary C.3. Fix s and t in [0, 1]. Then the mappings 

neN 

/rom W? )Odd (0, 1) x W^ even (0, 1) into W*(0, 1) and 

/ i — ^ ^^{cos[27ma; + 2s 2n (f)x] — cos 1-nnx) 

nGN 

/rom W / 2 s odd (0, 1) into W^O, -0 are analytic and Lipschitz continuous on bounded sub- 
sets. 

Indeed, for an odd / we have f{n) = —f(—n) and, as a result, 

s 2n (f) = | [/(-*) -/(*)] =i/(«); 

for an even g, we similarly have C2 n {g) = g(n) = g(— n )- Therefore the two series above 
coincide with ^$(2i/, g) + \g and the function ~/i(2i/) of Lemma ICT21 respectively, and 
the result follows. 
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Lemma C.4. For f and g in L 2 (0, 1), set 

*(/, 9) ■= V-P. E(- X ) n / 9{t) exp{[v™ + /(n)]i(l - 2t)} <ft e 2 ™*; 

then for every s G [0, 1] £/ie mapping 

* : L 2 (0,l)xW?(0,l)-). W 2 '(0,1) 
analytic and Lipschitz continuous on bounded subsets. 
Proof. The coefficient of e 27rma: in the above series for \l/ can be written as 

(C.l) f g(t) exp{/(n)i(l - 2t)}e" 2 ™^ 

Jo 

and gives the n-th Fourier coefficient of the function h := Y2T=ohk/k\, with ho '■= g, 
h%. := f <h> * M\g for h > 1, and Mi being the operator of multiplication by the 
function i(l — 2t), i.e., ty(f,g) = h. Since belongs to Wf(0, 1) and its norm there 
obeys the estimate 

ll^ll-^C^II/llo^llMillJllyll. 
with C being the constant of Proposition 113.21 and ||Mi|| s denoting the norm of the 
operator Mi in the space W 2 S (0, 1), we conclude that the mapping \I/ is analytic. Its 
Lipschitz continuity on bounded subsets is established in the usual manner. □ 

As above, by taking an odd / and an even g of zero mean, we see that the expressions 
of (IC.ll) for n and — n coincide with 

(C.2) 2(-ir / g(t)cos{[nn + s 2n (-if)](l-2t)}dt 



and give the 2n-th cosine Fourier coefficient of the function h =: *&(f, g) G W2 even (0, 1 
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